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O , Abstract 
<D 

c~| We study the toroidal dipole moment of the lightest neutralino 

in the constrained Minimal Supersymmetric Standard Model. The 

t—{ toroidal dipole moment is the only electromagnetic property of the 

neutralino. Since the neutralino is the LSP in many versions of the 
MSSM and therefore a candidate for dark matter, its characterization 
through its electromagnetic properties is important both for particle 
physics and for cosmology. We perform a scan in the parameter space 
of the cMSSM and find that the toroidal dipole moment is different 
from zero, albeit very small, in all the parameter space, and reaches 

H a value around 10~ 3 GeV~ 2 in a particular region of the parameter 

^s! space, well below experimental bounds. 

X 

a 1 Introduction 



One of the best motivated extensions of the Standard Model (SM) is the 
Minimal Supersymmetric Standard Model (MSSM), since, besides giving a 
solution to the hierarchy problem, provides us with a good candidate for cold 
dark matter (CDM), namely, the lightest neutralino. 

There are currently several experiments under way, and more planned for 
the future for direct and indirect detection of dark matter (DM) (for recent 
reviews on dark matter direction see flp]). If detected, it will be necessary to 
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discriminate between different candidates. To this end, it will be important 
to characterize as much as possible the different candidates. The neutralino 
is at present one of the best candidates for DM, and its electroweak properties 
can give us some insight into its nature. Because it is neutral, these properties 
appear only radiatively. 

Recently, there has been intense work on the electroweak properties of 
dark matter since they might be relevant in the calculation of DM decays 
and annihilations [3 15 , which have consequences in astrophysical processes 



16,17 and therefore are important in indirect astrophysical searches for DM, 



as in the calculation of the annihilation cross section of the DM itself. 

One of the least studied electromagnetic properties of a particle is the 
Toroidal Dipole Moment (TDM), which is directly related to the anapole mo- 
ment. The anapole moment corresponds to a T invariant interaction, which 



is C and P non-invariant 18 . It does not have a simple classical analogue, 



thus the toroidal dipole moment was introduced as a more convenient de- 



scription of T even, C and P odd, interactions 19 22 . The electromagnetic 



vertex of a particle can thus be expressed in a multipole parametrization, 
including the toroidal moments (see for instance |23|), which provides a one 
to one correspondence between the form factors and the multipole moments. 
Pospelov and ter Veldhuis have obtained an upper limit for the anapole 



moment of WIMPS 24 , using results from the DAMA and CDMS experi- 
ments 25 26 . In case the neutralino is the main component of dark matter, 
its anapole moment should comply with this limit. 

In this paper we calculate the TDM of the neutralino within the con- 
strained Minimal Supersymmetric Standard Model. We do a scan in the 
five parameter space of the cMSSM, and compare the results with the above 
mentioned experimental limit. The article is organised as follows: in sec- 
tion II we present a very brief summary of some aspects of the constrained 
MSSM (cMSSM) relevant to our calculation. In section III we review the 
general form for the electromagnetic vertex of a particle, and in particular 
for a Majorana particle, as the neutralino. We introduce the anapole moment 
and its relation to the toroidal dipole moment. In section IV we explain the 
methodology used to calculate the TMD of the neutralino in the cMSSM and 
we evaluate it for different values of the parameters. Section V presents the 
results obtained and our conclusions. 
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2 The MSSM and the neutralino as candi- 
date for dark matter 



The cold dark matter density is known to be 27 



n DM h 2 ~ 0.112, (1) 

where h is the Hubble constant in units of 100 km sec -1 Mpc -1 . The ther- 
mally averaged effective cross section times the relative speed of the dark 



matter particle, needed to get this relic density is 28 -30 

< ov >oc g A weak /lQ-K 2 m 2 x (2) 

consistent with the assumption of a weakly interacting dark matter particle 
(WIMP) with mass between 10 GeV - (few) TeV. 

The minimal supersymmetric extension of the Standard Model (MSSM) 
provides us with one of the best WIMP candidates for dark matter, namely 



the lighest neutralino (for reviews on SUSY see for instance 31,32]). The 
MSSM requires two complex Higgs electroweak doublets to give mass to the 
up and down type quarks and to avoid chiral anomalies. After electroweak 
symmetry breaking five physical Higgs states remain: two neutral CP invari- 
ant (h°,H°), two charged CP invariant (H + ,H~), and one neutral CP-odd 
(•>"). 

The MSSM has a new discrete symmetry, R parity, defined as R = 
(— l) 3 +2 + , where B and L are the baryonic and leptonic numbers re- 
spectively. This symmetry assignes a charge +1 to the SM particles and -1 
to the supersymmetric partners, thus making the lightest supersymmetric 
particle (LSP) stable. 

Supersymmetry has to be broken, or it would have already been observed. 
To break supersymmetry explicitly, without the reappearance of quadratic di- 
vergencies, a set of super-renormalizable terms are added to the Lagrangian, 
the so-called soft breaking terms. The Lagrangian for the soft breaking terms 
is given by 



Csoft = -^M a X a X a - -A ij %$ j( l) k - + c.c. - (m 2 )}0*Vi , (3) 

where M a are the gaugino masses, A^ k and B l i are trilinear and bilinear cou- 
plings, respectively, and (m 2 )* are scalar squared-mass terms. It is assumed 
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that supersymmetry breaking happens in a hidden sector, which communi- 
cates to the observable one only through gravitational interactions, and that 
the gauge interactions unifiy. This means that at the GUT scale the soft 
breaking terms are "universal", i.e., the gauginos M a have a common mass, 
as well as the scalars (m 2 )* and the trilinear couplings, A^ k . Requiring elec- 
troweak symmetry breaking fixes the value of and the absolute value of 
the Higgsino mixing parameter This is known as the constrained MSSM, 
or cMSSM, which is described by five parameters: the unified gaugino mass 
7711/2, the universal scalar mass mo, the value of the universal trilinear cou- 
pling A , the sign of Higgsino mass parameter /i, and the ratio of the vacuum 
expectation values of the two Higgses, tan (3. 

After the electroweak symmetry breaking the neutral and charged states 
in the MSSM can mix. In the case of the neutral ones they give rise to a set 
of four mass eigenstates, the neutralinos. It is the lightest one of these that 
is the LSP and a good candidate to dark matter in many SUSY models. The 
lightest neutralino, in the gauge eigenstate basis, is thus a function of the 
neutral Higgsinos and the neutral gauginos (Wino and Bino) 

iP = (B,W°,H u ,H° d ) . (4) 

The properties of the neutralinos will depend on the mixing, which in turn 
depends on the soft breaking parameters. Thus, the lightest neutralino can 
range from almost pure Bino to almost pure Higgsino. 



3 Toroidal Dipole Moment 

For 1/2-spin particles the most general expression for the electromagnetic 
vertex function, which characterizes the interaction between the particle and 
the electromagnetic field, is: 



(5) 



where /q(<? 2 ), f^q 2 ), fE^q 2 ) and /U(<? 2 ) are the so called charge, magnetic 
dipole, electric dipole and anapole form factors, respectively; where g M = 
v'u~Pn 1S ^e transferred 4-momentum; and a^ v = (i/2) [7^, 7, 
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These 

form factors are physical observables when q 2 — > 0, and their combinations 
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Figure 1: Current configuration with a toroidal dipole moment. The arrows 
on the torus indicate the direction of the current, and the TDM is directed 
towards the symmetry axis of the torus. 



define the well known magnetic dipole (/i), electric dipole (d) and anapole 
(a) moments. 

The electromagnetic properties of Majorana fermions (like the neutralino) 
are described by a unique form factor, the anapole, /a(<? 2 )- This is a con- 
sequence of CPT-invariance and the C, P, T properties of T At (g 2 ) and the 
interaction Hamiltonian. Thus, the electromagnetic vertex function of a neu- 
tralino can be writen as 

r,(q 2 ) = f A (q 2 )(q 2 l,-^H- (6) 

The anapole moment was introduced by Zel'dovich to describe a T- 
invariant interaction that does not conserve P and C parity [18] . The anapole 
moment does not have a simple classical analogue, since Ja{q 2 ) does not cor- 
respond to a multipolar distribution. A more convenient quantity to describe 



this interaction was proposed by V. M. Dubovik and A. A. Cheshkov 34 
the toroidal dipole moment (TDM), r(q 2 ). 

The TDM and the anapole moment coincide in the case of = rrif, i.e. 
the incoming and outgoing particle are the same. This type of static multi- 
pole moments does not produce any external fields in vacuum but generate a 



free-field (gauge invariant) potential 21 , which is responsible for topological 
effects like the Aharonov-Bohm one. 

The simplest TDM model (anapole) was given by Zel'dovich as a con- 
ventional solenoid rolled up in a torus and with only one poloidal current, 
see fig. [TJ For such stationary solenoid, without azimuthal components for 
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the current or the electric field, there is only one magnetic azimuthal field 
different from zero inside the torus. 

4 One-loop calculation 

The TDM of the neutralino may be defined in the one-loop approximation 
in the cMSSM by the Feynman diagrams shown in figs. [2] and |3j where 
/ represents the charged fermions of the SM. Taking each fermionic family 
separately we obtain 94 Feynman diagrams in total: 66 corresponding to 
self-energy and 28 to vertex corrections. 




Figure 2: One-loop vertex corrections to the process 7 — > XiXv 

We use FeynCalc to calculate the amplitude of these diagrams. Since we 
are only interested in the terms that contribute to the anapole form factor, 
we isolate the ones that have the Lorentz structure 7^75. It is important to 
notice here that we work in the t'Hooft-Feynman gauge (£ = 1). One of the 
first results we obtain is that the self-energies 'jH , 7/1 , jA and 7G do 
not contribute to the calculation. If we call Sj the coefficient that multiplies 
7ju 75 f° r the ith diagram, then we have that 

5> = /a(?V- (7) 
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Figure 3: One- loop corrections to the self-energy for the process 7 — > XiXv 



To obtain the toroidal dipole moment r = /a(0) we use the l'Hopital rule 
and get 

r = f A (Q) = lim = | 92 _ . (8) 

g 2 ->o q 2 oq 2 

The contributions to the self-energies have two point Passarino-Veltman 
scalar functions of the type B (q 2 ,x 2 ,x 2 ) and B (0, x 2 , x 2 ). Likewise, the 
contributions to the vertex corrections have two and three point scalar func- 
tions of the type B (g 2 , x 2 , x 2 ), B y 2 , a; 2 ) and C (q 2 , M| ? , M| ? , x 2 , x 2 , y 2 ) 

In both and y represent the masses of the particles in the loop. 

When evaluating (6), derivatives of the Passarino-Veltman functions ap- 



pear. To evaluate the -Bo's, as well as their derivatives, we use LoopTools 35 
To evaluate the Co's and their derivatives we expande them in a power series 
around q 2 = (see appendix). 

The expression obtained for the toroidal dipole moment depends on vari- 
ous parameters of the MSSM, including the supersymmetric particles masses 
as well as the mass mixing matrix elements, the value of tan /3, and the values 
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of the soft breaking terms. We evaluate the TDM within the cMSSM using 
Suspect [36 1, by fixing the value of A , tan/3 and sign//, and scanning over 
the other two parameters mo and mi/2, from to 1500 GeV and 250 to 1500 
GeV, respectively. 

Figure [4] shows the neutralino toroidal dipole moment for tan/3 = 10, 
fi > 0, and three different values of A , —1000, and 1000 GeV (top to 
bottom). Comparing the three different plots, no dependence on Aq is shown. 
The TDM is very low for almost every region of the parameter space scanned, 
with values between 10~ 5 and 10~ 8 GeV -2 . However the TDM increases for 
increasing tuq and decreasing m^, reaching values over 10~ 3 GeV -2 for high 
m o (> 800 GeV) and low (< 400 GeV). Similarly, figure [6] shows the 
results for tan/3 = 50, fi > 0, and the same three different values of A , 
— 1000, and 1000 GeV (top to bottom). The TDM reaches values around 
10- 3 GeV~ 2 . 

Figure [5] shows the TDM for tan j3 = 10 and /i < 0. Sign/i > may 
solve the problem of the discrepancy between the measured value of g — 2 
of the muon and the one predicted by the SM. However, this does not mean 
negative sign/i is ruled out since others mechanisms could solve this problem, 
and therefore sign/i < should be taken into consideration. 

Figure [7] shows a comparison of two plots for different tan but same 
sign/i and A . This figure shows the dependence of the TDM on tan/3. 
Figure [8] shows a comparison of two plots for different values of sign/t but 
same tan/3 and Aq. This figure shows no dependence of the TDM on sign/i. 

Notice that in all the plots the region for which M^o = Mf is suppressed 
since we are not considering this possibility. This condition (M^o = M f ) 
separates the region where the neutralino x? is the LSP and the one where 
the stau f is the LSP. 

If combined with other criteria, our results can be used to further reduce 
the parameter space of the cMSSM. As an example, figure [9] combines our 
results with those obtained in ref. [37]. They studied the allowed parameter 
space for diferent models (including cMSSM), supposing that the neutralino 
is the main component of the CDM, varying Aq and imposing the restrictions 
of cosmological experimental data from WMAP. As can be seen in the plot, 
their criteria favours a region with very low m , which is compatible with a 
TDM < 10~ 5 GeV -2 for the neutralino. In other words, if a TDM is measured 
for a WIMP and it is around or higher than 10~ 4 GeV, the neutralino of the 
cMSSM would no longer be a good candidate for CDM in this region of 
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parameter space. 

In figure 10 we compare our results with a frequentist analysis of the 
probable ranges of parameters of the MSSM using the results of initial di- 
rect LHC searches for supersymmetry, combined with the required cold dark 
matter density, and the spin-independent dark matter scattering cross sec- 
tion [38]. The figure shows the 68 and 95% CL contours (red and blue, 
respectively) both after applying the CMS and ATLAS constrains (dashed 
and solid lines, respectively) and beforehand (dotted lines), taken from their 
paper. Also shown in the figure as open (solid) green stars are the best-fit 
points found after applying the CMS (ATLAS) constrainsts in the cMSSM. 
As can be seen in the figure, this analysis also favours a neutralino with a 
very low TDM. 

In figure 11 we compare our results with another frequentist analysis 
of the CMSSM pararameter space including the public results of searches 
for supersymmetric signals using 1/fb of LHC data recorded by ATLAS 
and CMS and 0.3/fb of data recorded by LHCb in addition to electroweak 
precision and B-physics obsevables [39]. The figure shows the 68 and 95% 
CL contours (red and blue, respectively) with LHCi/f, data (solid lines) and 
showing pre-LHC fits (dotted lines), taken from their paper. Although this 
analysis widens the allowed parameter space, it still favours a neutralino with 
very low TDM. 



5 Conclusions 

We calculated the only electromagnetic property of the lightest neutralino: 
its toroidal dipole moment. Its characterization is extremely valuable for dis- 
criminating different models which have the neutralino as dark matter candi- 
date. We performed the calculation in the framework of the cMSSM, however 
a similar analysis can be performed for other models (work in progress). We 
found that the TDM of the neutralino is highly sensitive to m , mi/2 and 
tan (3, but very weakly or practically non-dependent on A and sign//. 

All points in the parameter space we scanned give a TDM consistent with 



the upper limit (~ 1CT 2 GeV -2 ) obtained by Pospelov and ter Veldhuis 24 
for WIMPs interacting with heavy nuclei using data from the CDMS and 
DAMA experiments. However, this data can and will be improved in the 
next few years helping to refine the upper limit, likely ruling out some regions 
of the parameter space. 
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The TDM analysis can be used as another criteria to constrain the pa- 
rameter space of a given model which has a neutralino as candidate for dark 
matter. Thus, according to our results, if a non-zero (around 10 _4 -10 -3 
GeV -2 ) TDM could be measured for the neutralino, that would indicate 
that the favored region of the parameter space of the cMSSM would be high 
m o (> 800 GeV) and low mw2 (< 400 GeV). Otherwise, other regions are 
compatible with a TDM lower than 10 -5 GeV. 

If combined with other criteria (such as dark matter relic density, or fre- 
quent ist analysis of data from the LHC) the parameter space can be reduced 
even further. In fact, the combination with these other criteria favours a 
region with low m , which is compatible with a 10~ 5 -10 -8 GeV -2 TDM for 
the neutralino. This means, among other things, that if a TDM higher than 
10 -4 GeV~ 2 is measured for a WIMP it would exclude some regions of pa- 
rameter space of the cMSSM (and other more specific models), at least if the 
neutralino is the only component of dark matter. One should keep in mind 
that all analysis mentioned were made with this assumption. 
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Appendix: Scalar Three-point Function 

In this appendix, we analyse the Passarino-Veltman scalar three-point func- 



tion Co(q 2 , x 2 , x 2 , z 2 , z 2 , y 2 ) (40 41 which appears in the TDM calculation. 
Here q 2 denotes the photon transfered 4-momentum, x is the neutralino mass, 
and y and z are the masses of the particles running in the loop. 



The corresponding plot for this Co function can be seen in figure 12 
The red line shows the numerical solution, the blue line line represents the 
approximate solution, i.e., the Taylor expansion around q 2 = 0, which can 
be written as follows: 

C (g 2 , x 2 , x 2 , z 2 , z 2 , y 2 ) = a + a x q 2 + 0(g 4 ). (9) 
The coefficients a, are functions of the masses: 
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log 

a ° = — 2cc 2 / +aiogu > ( 10 ) 



x 4 — y 2 x 2 — 2z 2 x 2 + z A — y 2 z 2 1°S v , ^ 

1 6x 2 z 2 (—x + y — z)(x + y — z)(—x + y + z)(x + y + z) 12x A ^ ' 

(11) 

where 

[ix 2 + iy 2 - iz 2 + -J -y 4 + 2(a; 2 + z 2 )y 2 - (z 2 - a; 2 )) (ia; 2 - iy 2 + iz 2 + -J -y 4 + 2(ai 2 + z 2 )y 2 - (z 2 - x 2 )j 

(~ix 2 + iy 2 - iz 2 + ^/-y 4 + 2(a; 2 + z 2 )y 2 - (z 2 - a; 2 )) (-ix 2 - iy 2 + iz 2 + -y 4 + 2(x 2 + z 2 )y 2 - (z 2 - a; 2 )) 

(12) 



i(x 2 + y 2 — z 2 ) 



2x 2 y/~- 



2y 2 x 2 + 2z 2 x 2 -y A -z A + 2y 2 z 2 



and 



(13) 



i{ x 2 + y 2 - z 2 )(x 4 



■4y 2 x 2 -2z 2 x 2 + y 4 + 4 -2y 2 z 2 ) 



I2x 4 (-x + y- z)(x + y- z)(-x + y + z){x + y + z) y '-x 4 + 2y 2 x 2 + 2z 2 x 2 - y 4 - z 4 + 2y 2 z 2 



(14) 
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Figure 4: TDM for tan/3 = 10, sign/i = + y A = -1000 (top), (centre) y 
1000 (bottom) GeV. 






Figure 5: TDM for tan/3 = 10, sign/i = — y A = —1000 (top), (centre) y 
1000 (bottom) GeV. 






Figure 6: TDM for tan/3 = 50, sign/i = + y A = —1000 (top), (centre) y 
1000 (bottom) GeV. 
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Figure 7: TDM for sign/i = +, A = GeV, tan/3 = 10 (green) and tan/3 = 
50 (red). 



sign(^i)=+1 + 
sign(u)=-1 




Figure 8: TDM for tan /3 = 10, A = GeV, sign/i = + (red) and sign/i = — 
(green) . 
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Figure 9: TDM results for tan = 10, Ao = GeV and sign/i = + combined 
with restrictions from cosmological observational data (WMAP). 
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Figure 10: TDM results for tan /3 = 10, Aq = GeV, sign/x = + combined 
with restrictions from initial direct LHC searches for supersymmetry. 
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Figure 11: TDM results for tan/3 = 10, Aq = GeV, sign/x = + combined 
with restrictions from searches for supersymmetric signals using 1/fb of LHC 
data. 
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Figure 12: Comparison between numerical (red line) and approximate (blue 
line) scalar three-point function Co(q 2 , x 2 , x 2 , z 2 , z 2 , y 2 ), with x = 97.7 GeV, 
y = 415.4 GeV and z = 80.43 GeV. The analytical approximation (blue line) 
is only valid for q 2 — > 0. 
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